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Abstract. Originally introduced in solid state physics to model amorphous materials and 
alloys exhibiting disorder induced metal-insulator transitions, the Anderson model — 
—A + VL on / 2 (Z d ) has become in mathematical physics as well as in probability theory a 
paradigmatic example for the relevance of disorder effects. Here A is the discrete Laplacian 
and K, = {Vu,{x) : x G Z d } is an i.i.d. random field taking values in R. 

A popular model in probability theory is the parabolic Anderson model (PAM), i.e. the 
discrete diffusion equation dtu(x,t) = —H LO u(x,t) on Z d x R+, u(x,0) = 1, where random 
sources and sinks are modelled by the Anderson Hamiltonian. A characteristic property of 
the solutions of (PAM) is the occurrence of intermittency peaks in the large time limit. These 
intermittency peaks determine the thermodynamic observables extensively studied in the prob- 
abilistic literature using path integral methods and the theory of large deviations. 

The rigorous study of the relation between the probabilistic approach to the parabolic An- 
derson model and the spectral theory of Anderson localization is at least mathematically less 
developed. We see our publication as a step in this direction. In particular we will prove 
an unified approach to the transition of the statistical moments (u(0, t)) and the integrated 
density of states from classical to quantum regime using an effective medium approach. As a 
by-product we will obtain a logarithmic correction in the traditional Lifshitz tail setting when 
Vui satisfies a fat tail condition. 



1. Introduction 

The Anderson model is the family of discrete random Schrodinger operators {H^} defined by 

H u = -A + V u . 
Here A is the discrete Laplacian on Z 2 (Z d ) 

[A«](x) = Yl l u (v) ~ u ( x )]- 

\x-y\=l 

The random potential {Vu}(x)} x£Z d is a field of independent and identically distributed random 
variables with common distribution Pq. Denoting the expectation value by ( . ) we assume 

G(t) := log(exp(-tK;(0))> < oo (1) 
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for all t > . {H w } is an ergodic family of self adjoint operators on l 2 (Z, d ). In many con- 
crete situations exponential localization is proven at the bottom of the spectrum [15, 26, 30], i.e. 

• dense point spectrum close to inf a{H u ), 

• exponentially decaying eigenfunctions. 

The spectral analysis of {H^} is motivated by applications in solid state physics, e.g. local- 
ization phenomena, electrical resistance, low temperature physics, We refer to [19] and 

references therein. 

The parabolic Anderson model (PAM) is the discrete diffusion equation with random sources 
and sinks: 

d t u(x,t) = -H w u(x,t) (x,t) £Z d x [0,oo), 

u(x,0) = 1 x£Z d . 

Assuming (1) the parabolic Anderson model has a.s. an unique, nonnegative solution given by 
the Feynman-Kac-representation [7] 



u(x,t) = E x 



(2) 



Here ¥?[.] is the expectation value of the random walk in continuous time generated by —A 
starting in x. For t > the random field {u(x, t) : x £ is stationary, ergodic and mixing 
under translations. The moments (u(0,t) p ) and the correlation function are finite [7, 12]. 
Describing the large time diffusive behaviour of a classical particle in a random medium with 
traps the applications of the parabolic Anderson model are numerous. (PAM) is used as a 
linearised model of chemical reaction kinetics exhibiting macroscopic pattern formation in the 
spatial distribution of reagents, has interpretations in polymer physics and is used to describe 
population dynamics in an inhomogeneous environment modelling the availability of nutrients. 
For a very recent application of (PAM) in this biological setting as well as for a comprehensive 
summary of other interpretations, respectively interesting generalizations of (PAM) we refer 
to [20] and references therein, see also [12, 7, 10, 25]. 

In the limit t — > oo the solution u(x, t) shows a.s. a very strong spatial inhomogenity caused 
by very rare potential constellations. This phenomen is known in the probabilistic literature 
as intermittency and is described by asymptotic behaviour of the moments {u(0,t) p ) [7]. As- 
suming V UJ (x) > the first moment {u(0,t)) can be interpreted as the survival probability of 
a particle that is put randomly on 7L d . 

The intuitive link between the probabilistic and the spectral point of view is: 

Shape of intermittency peaks < — > Localized eigenfunctions, 
Local killing rate < — > Eigenvalues. 



A quantity to formalize the intuitive link between the Anderson model and PAM is the the 
integrated density of states measure v ( [15, 17, 31] and references therein). Here we are 
interested in the integrated density of states (IDS) N(E), i.e. the distribution function of v 

N(E) := u((-oo,E]) = lim |A| "^{eigenvalues of H% < E} (3) 

|A|— s>oo 
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with 



XaHujXA- 



(4) 



The integrated density of states N(E) is the fundamental quantity to study the thermodynam- 
ical properties of disordered systems. Moreover, N(E) is used to prove localization properties 
of the system. In particular we are interested in Lifshitz tails, i.e. the behaviour of the IDS in 
the limit E \ inf (j(H u ). Assuming (1) the Laplace transform 



N{t) 



-xt 



dv{\) < oo (t > 0) 



of v exists [14] and has the Feynman-Kac representation ([4], [14]) 



N(t) = (E° 



exp 



V UJ {x s )ds 5 (x t ) 



(5) 



(6) 



The first proof of Lifshitz behavior (for the Poisson model ) was given by Donsker and Varadhan 
[6]. Starting from the Feynman-Kac representation their estimate of N(t) in the limit t — > oo 
relied on an investigation of the "Wiener sausage" and the machinery of large deviations for 
Markov processes developed by these authors. To obtain information about the behavior of 
N(E) for E \ inf ^(i?^) from the large t behavior of N(t) one uses Tauberian theorems [3], 
see also Appendix 2. This technique was already used by Pastur [1, 27]. The behaviour of 
N(t) in the limit t — > oo. is also closely related to the long time behaviour of the moments 
(it(t, 0)) of the parabolic Anderson model. 



To formulate our main result Theorem 2, we remind the definition of regularly varying functions 
and of the de Haan class [3], see also Appendix 1. 

Definition 1. 

(i) A function g > defined on some neighbourhood [X, oo) of infinity satisfying 

g{\t)/g(t) ^°A"(l + o(l)) 

for all A > is called regularly varying of index p. We write g G R p . If p = then g 
is said to be slowly varying. If g varies regularly with index p, we have g(t) = t p go(t), 

(ii) For g G R p and A £ (0, 1] the de Haan class U g is the class of functions H : R — > M 
satisfying 

H{t) - H(Xt) c g h p {\)g(t){l + o(l)), 
where g G R p is called the auxilary function and c g is the g-index. 

Our main result estimates the Laplace transform N{t) defined in (5) and the first moment 
{u(t, 0)) in terms of two variational functionals. Here u(t, 0) is the solution of the parabolic 
Anderson model. The variational functional of the lower bound is given by 

X J(t) := Adsin 2 ("j^y) + c g h p (l~ d )g(t) i G N (7) 
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and 

r max 1/2 <h<i min [2d(l - 2y/l^h)^/2 + (1 - h)c g h p (l - h)g(t)} 1 = 1, 

Xt(t) ■= I o . i v (8) 
7 sin 2 (f jljj +4' 1 c 5 / lp (4£)- d ) 5 (t) £ > 1, 

£ G N, 7 > 0, is the corresponding variational functional of the upper bound. 

Theorem 2. Suppose G(t) < oo, t > and G(t)/t G IT 9 lozi/i auxilary function g(t) G .Rp, 
p G [— 1, oo), g-index c g and tg(t) — > oo m t/ie Zimif i — >■ oo. Denote by N(t) the Laplace 
transform of the integrated density of states and by (u(t, 0)) the first moment of the solution 
of the parabolic Anderson model. Then with xf(t) as defined above 

G(t) - t inf X 7(i)(l + o(l)) < \ogN{t) < log(u(t, 0)) < G(t) - t inf x +(t)(l + o(l)). (9) 

Remark 3. Due to 0)n(s, 0)) = + s,0)) [7] Theorem 2 can also be used to estimate 
the higher moments of u(t, 0). 

Theorem 2 is motivated by the theory of critical phenomena in statistical physics. Looked at 
from this angle the variational problem in (9) correponds to the minimization of the free energy. 
The ground state energy of the Dirichlet-Laplacian on = A^(0) := {x G 7L d : |x|oo < £} given 
by 

'n 1 
sin" | — - 



2£+ 1 



plays the role of an order parameter. The parameter I corresponds to the extension of the 
Lifshitz ground state. The dependance on the random potential is encoded in the effective 
potential G(Xt)/Xt, X G [0,1], respectively the deviation 

S(X,t) := G(t)/t - G(Xt)/Xt > (10) 

to the maximal effective potential value G(t)/t. Denoting by £*(t) the optimal length defined 
by 

we can distinguish two qualitatively different regimes in dependence on the single site distri- 
bution [26]. 

Quantum regime: If we assume that G(t)/t G II g with 

S(i)^°0 (11) 

the asymptotics of (u(t,0)) and N(t) are dominated by the energy form of the Laplacian, i.e. 
by the quantum kinetic energy, respectively the rate function of the occupation time measure 
in large deviation theory. As a consequence £*(t) will tend to infinity and from the physical 
point of view Lifshitz tails are expected [22], i.e. an asymptotic behaviour of the IDS like 

N(E) ~ C ie - C ^ E - E ^~ d/2 . (12) 
The content of the next Corollary is that this is only approximately correct. 
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Corollary 4. Suppose G(t) < oo, t > and G(t)/t G II 9 with auxilary function g(t) = 
t p 9o(t) ^ Rp, P ^ [ — 1, 0] , g-index c g and go G Rq s.t. 

9o(tgo(t) 1/p )/9o(t) ^° 1. (13) 

Furthermore assume that g(t) — > and i<?(i) — >■ oo in the limit t — > oo. T/ien £/ie asymptotic 
optimal length is given by 

t(t) g(t)V(*»-2) 

and t/iere exists constants C±,C2 > 0, s.f 

- c x tt(t)- 2 < logiv(t) < iog(«(t,o)) < -c 2 tt(t)- 2 . (14) 

Suppose that p G [—1,0). TTien i/ie integrated density of states satisfies 



-C 1 E-I +1+ ' , ~ 1 g (CE^-^P) 1/P < log N(E) < -C 2 E~i +1+ ^ gJcE^^^ 



(15) 

with C,Ci,C 2 >0 and E\0. 

Choosing for example the uniform distribution in [0, 1] we have 

M _ ffifi 6 ^ (i 6 ) 

i At At K J 

and (15) becomes 

CiE~ hog E < log N(E) <C 2 E~ hog E (17) 
in the limit E \, 0. Comparing (17) and the estimate 

log(-log(iV(£))) = _d 
E\o log(-E) 2 

proven with spectral theorectic methods in [16, 29] assuming the fat tail condition 

P[K,(0) < -E] C£ fc , A; G N , (18) 

we obtain a logarithmic correction predicted in the physics literature [23, 28]. 
Assuming p G [1,0) the assumptions of Corollary 4 corresponds in the probabilistic setting of 
[4] to the existence of a non-decreasing function t i— >■ at G (0, oo) and a function G: [0, oo) — )• 
(— oo, 0], G j£ 0, such that 

a d+2 / t \ 

li m _^ G y )=G(y), y > 0, (19) 

uniformly on compact sets in (0,oo) (p = is discused in [11]). Condition (19) is satisfied if 

F[V u (0)<E] exp(-C J E" £ 7 i ). 

Using the Feynman-Kac-representaion (2) and the large deviation theory for path integrals it 
is proven in [4] that 

i logMO, tf) ^ - \ ( X + o(l)) , (20) 

Pt pta/ a z pt 
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with 

X 



inf {I iVgix^dx-Cp' 1 [ g{x) 2{1+p) - g(xf dx), (21) 



II S II 2 = 1 



An application of Tauber theory gives in the limit E \ 

logiV(E) = C(p,x)E~-^ +0 ^. 

Finally let us discuss the almost bounded single site distributions, i.e. G[t)/t £ II 3 with g £ i?o, 
lim^oo gf(t) = 0. This setting is again dominated by the kinetic energy and £*(t) — > oo. The 
probabilistic counterpart of (14) 

1 G(pta~f) 1 . In 

- log(u(i, Of) = {P f/ - — (pd(l - - log P.) + oil)) , (22) 
pt pta pt d a z pt 2 vr 

as t — > oo is proven in [11]. Here the scale function at is defined by the fixed point equation 

g{taT d ) = a~ t 2 . (23) 

Furthermore we want to refer to [18], where in Theorem 1.5 for unbounded single site distri- 
butions satisfying G(t)/t £ H g with g £ Rq, lim^oo g(t) = the asymptotics of the IDS in the 
limit E — > — oo is proven. 

Classical regime: Let us now consider the classical regime, i.e. G{t)/t £ Ii g with 

liminf g(t) > 0. (24) 

t— >oo 

Then the quantum kinetic energy/occupation time measure and the effective medium are on 
the same scale, respectively g(t) dominates the energy form of the Laplacian. As a consequence 
£*(t) stays finite in the limit t — > oo and the IDS is given by the shifted rate function of the 
single site potential. Let us first discuss the asymptotics of the statistical moments and of 
N(t). 

Corollary 5. Suppose G(t) < oo, t > 0, G(t)/t £ U g with auxilary function g £ R p , p £ [0, oo) 
and g-index c g , s.t. (24) is satisfied. With 



X-{t): 

and 

X\{t) : 



fl c g g(t) >2vr 2 

\4cgg(t) + c g g(t) log (2vr 2 /{c g g(t))) c g g{t) < 2tt 



' 1 - 2{c g g{t)r 1 / 2 c g g{t) > 2e 2d + 7 n 2 /2d, 

min [ 7 /(4d), + ^ log (gg)] c g g(t) < 2e 2d + 7^2/2^. 

we have 

G(t) - 2dt X *-(t){\ + o(l)) < log < log(u(t, 0)) < G(t) - 2dt X +(t)(l + o(l)). (25) 

While g(t) dominates the variational functional for p > 0, the slowly varying functions define 
the borderline between the quantum and the classical regime. The main objective of [8] is the 
double exponential distribution 

P[K,(0) < E] E ^J°° exp (-e^M (26) 
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with G(t) = c g t log(c 9 t) — Cgt + o(t), respectively S(X, t) = c g log(A) + o(l). The probabilistic 
approach obtain (20) with 

X= min {- V (g(x) - g(y)) 2 - p V g 2 (x) log g 2 (x)\ . (27) 

\x-y\ = \ 

and ~ l/y/c g G (0, oo), i.e. the intermittency peaks are finite but nontrivial. As a conse- 
quence there are no Lifshitz tails. The IDS is the single site rate function 

I(E):=M[Et + G(t)\ (28) 

shifted by the constants 2d\± encoding the size of the intermittency peaks. 

Corollary 6. Suppose G(t) = c g tlog(c g t) — c g t + o{t) and x± Q£ in Corollary 5. Then 

-Cl{E-2d X *.)(l + o(l)) < log N(E) < -l(E-2d X *+)(l + o(l)). (29) 

Finally let us discuss single site distributions satisfying G(t)/t G II 5 with g G R p , p > 0, i.e. 
lim^oo g(t) = oo. An example is the Weibull distribution 

F[VU0)<E] ex V (-C(-E) a ), 

a > 1. The asymptotics of (u(t,0)} and N(t) are dominated by g(t). As a consequence we 
have £*(t) = 1 and the asymptotic behaviour of the IDS is given by the maximal shifted single 
site rate function. 

Corollary 7. Suppose G(t)/t G II S with g G R p , p > 0. Then 

log N(E) = -I{E-2d){l + o{\)). (30) 
Corollary 7 corresponds to the results obtained in [9, 18]. 

To prepare the discussion of our strategy to prove Theorem 2 let us mention the key ideas of 
the spectral and the probabilistic argumentation. 

Optimal-Fluctuation Method The core of the spectral theoretic approach is a rare region effect 
[32] predicted by Lifshitz [21, 22]. Let us assume that > and inf(r(ff w ) = 0. To find an 
eigenvalue smaller than E > 0, the uncertainty principle forces the potential V u to be smaller 
than E on a large set whose volume is of order E~ d / 2 . This is a very rare event and its 
probability is approximately 

P[Jj{x G A : V{x) <E}> E~ d ' 2 } w e -CE~ d /\ (31) 

Applying Dirichlet-Neumann-bracketing, perturbation theory or periodic approximation the 
heuristic argument above can be proven rigorously. 

Path-integral and large deviation techniques The probabilistic methods combine the Feynman- 
Kac-representation of (u(0,t)) and N(t) and large deviations techniques for path integrals. 
Informally the key idea is to express (u(t, 0)) by means of local times of random walks on Z d 

k(x) := [ l Xs=x ds xeZ d ,t>0, (32) 
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respectively the occupation time measure Lt := lt{x)/t and to average w.r.t. the random 
potential 



(u(x,t)) = (K x [exp(J2 V u (x)k 



x£Z d 



E x [expG(t) + t ^ \[G{L t {x)t) - L t {x)G{t)})]. 



t 

The next step is to represent the expectation value above as a Laplace integral for the occu- 
pation time measure, to apply large deviation principles and Varadhan's Lemma to obtain 

(u(x,t)) 

= exp(G(t) + o(t)) f e~ l ^ d ^G(L t (x)t)-L t (x)G(t)] £ ^ 

= e W (G(t)-ta(t)- 2 X (l + o(l))) 

as t — > oo. Nevertheless a rigorous implementation of the argument above is nontrivial (a good 
guess of at is necessary) and has to be proven in four steps: 

• making the space finite (but still time-dependent), 

• using a Fourier expansion and scaling properties, 

• removing the time-dependence of the box (compactification) , 

• applying the large deviation arguments. 

As a consequence of the technical problems resulting from the mathematical implementation 
at least up to now an unified approach treating all single site distributions at once does not 
exist. 

The strategy used in the proof below combines ideas from spectral and probability theory. The 
toehold proven in Section 4 is to restrict the estimates of {u(x, t)) and N(t) to a cube A = A(i) 
of time-dependent side length L = L(t) and to study 

(exp (-t E x (Ha))) = ( sup exp(-t [(^|-AiVp) + (v^|WP)])). (33) 

peart! (A) 

with Ei(H^) = inf a(H®) and the set of probability measures with support contained in A 

arti(A) := {p6 £W 1 (Z d ) : supp^C A}. (34) 
In (33) two competing effects are coupled: 

• High productivity in (■s/p\V u} ^ s /p) with respect to V u 

versus 

Small probability of extreme productive values of the potential 

• High productivity in (■ s /p\V u} ^ s /p) with respect to p 

versus 

Small probability of the occupation time encoded in (y/p\ — A^^/p). 
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To prove Theorem 2 we want to use, that the optimal p balancing between the two competing 
effects above satisfies 

• p is concentrated on a small cube C A, 

• p is relatively uniform on A^. 

The proof of the lower bound of Theorem 2 in Section 2 is elementary. We can interchange 
the supremum and the expectation value in (33) and obtain an effective medium problem. 
By restricting to a subset D C 9J?i(A) and optimizing with respect to 53 we obtain the lower 
bound. 

The upper bound of (33) proven in Section 3 is slightly more difficult. We have to control all 
p £ 9JTi(A) and it is not possible to interchange the supremum and the expectation value in 
(33). The first step in the proof is the classification of p £ SOti(A) in Definition 13 below. From 
the spectral theoretic point of view Definition 13 corresponds to a classification with respect to 
the kinetic energy while from the stochastic point of view the classification is a down to earth 
variant of the contraction principle concerning the asymptotic probability of the occupation 
time measure p. The problem is then to estimate 

( sup exp ( ty2p(x)V(x) j ) (35) 

solved by an effective medium theory. Finally in Section 4 we prove that all occuring error 
terms are negligible compared to the first correction of G(t) given by tinf^g^x^^)- 
Let us summarize the current state of research. We discussed quite at lot of publications 
based on probabilistic and on spectral theoretic methods. The spectral approach is close to 
the physical intuition, but the assumptions are restrictive. Moreover important aspects of the 
phenomenology get lost. This concerns the dependence of the IDS on the single site distribu- 
tion, e.g. the logarithmic correction for fat tail distributions. The probabilistic approach deals 
all single site distributions and obtain sharp asymptotics. Nevertheless an unified approach 
systematically explaining the relevant effects like the transition from the quantum to the clas- 
sical regime does not exist. Symptomatically the probabilistic publications are motivated by 
spectral theorectic heuristics, while the goal of the formal proof consists in guessing a good 
scale function, s.t. a large deviation principle can be applied. The intrinsic motivation of the 
scaling remains unclear. 

While not obtaining the sharp asymptotics at least partially Theorem 2 resolves some of the 
questions discussed above. The problem is discussed in an unified setting. The distinction 
between quantum and classical regime is a natural consequence of the variational description. 
The elaborated large deviation techniques for path integral measures as well as the preparations 
to apply them are avoided. Finally an important motivation for our approach is understanding 
the structural relation between the probabilistic and spectral theoretic methods. 

2. The lower bound 

Assuming the hypotheses of Theorem 2 we want to prove in this section the following lower 
bound. 

Proposition 8. With S(.,.) as in (10) we define 

X7(t)=4dsin 2 (lJ-^ + s(r d ,t). 
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Then 



(u(t,0)) > N(t)> explG(t)-tinfx7(t)](l + (l)). 

Lemma 9 below is strongly influenced by the probabilistic strategy of commutating the expec- 
tation values of the random potential and of the random walk to obtain an effective description. 

Lemma 9. Denote by E\(H®) = inf a(H^) the ground state energy of defined in (4), 
A = Aj(0) := {x G Z d : \x\oo < I}. Then 

(exp (-i > exp (c(t) - ijmtf ^ ^ (y/p\ - Af ^) + J p(x)5(p(x), t)^ . 

Proo/. 

<exp(-t E 1 (#£))> 

> sup <exp(-t (vp|-Afvp) + (v^lK,VP))) 
peOTi(A) 

= sup exp(-t - A^^/p)) TT(exp(-tp(i)K ; (ac))) 

P ewt 1 (A) xeA 

= sup exp I -i (vp| - AfvP) + J2 G(p(x)t) ) 
P ewi(A) V xGA / 

= exp^-t^inf^^Vpl-AfvP) + &W%),t)|J. 



□ 

The next problem is to distribute the probability mass = 1 of the occupation time measure 
p £ 9Jli(A) s.t. the competition between diffusion and particle creation encoded in 

(Vp| - Af y/p) + ^p(x)5(p(x),t) (36) 
xeA 

is minimized. We have to balance between: 

• The energy form of the discrete Laplacian 

(VpI&aVp), (37) 

i.e. the rate function of the occupation time measure encoded in p G 9JTi(A) [12]. If 
t > is large, it is much more likely that the local time is smeared over a large region 
than being localized in small subset of A. 

• The second term in (36) expresses the opposite effect. Due to convexity of G(t) we 
have 

EW) ^ = G(t), (38) 

xGA xGA 
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respectively 

Y, P {x)S(p{x),t) > (39) 

x£A 

for general p € 9Jti(A) and equality if there is a x 6 A s.t. p = 5 X . The function S(X,t) 
measures the deviation of the productivity of p £ 2Jti(A) compared to the maximal 
productivity given by G(t). 

To deal the competition between diffusion and particle creation s.t. the lower and upper bound 
are in good agreement we restrict 9Jti(A) to the subset 2) C £PTi(A) defined below. 

Definition 10. Denote by — , 1 < I < I, the restriction of the discrete Laplacian to the 



box A e = A e (0) := {x £ Z d : \x\oo < £}, by 

4>l : A e -> [0, oo) 



its ground state and by 



*W = n(7TT)*-(^) 

J'=l 



its ground state energy. Then, 

D := {<$ : l<£< I}. 

The set D contains relatively uniformly distributed prototypes of occupation time measures 
localized in a small volume. Inserting these intermittency peak candidates in Lemma 9 we 
obtain the following estimate. 

Proposition 11. Denote by E\(H?) = inf o~(HP) the ground state energy of . Then, 

(exp (-t E 1 {Hi))) > exp (o(t) -tM^xT(t)\ . 

Proof. Lemma 9 and Definition 10 yield that 
<exp(-i E 1 (iJf))) 



> 



exp \G{t) -t inf y (Vp\&aVp) + Y f p(x)S{p(x),t) 



In particular choosing yjp = <f>i with <j>g : A^ — > [0, oo) and E\ y— A^J as in Definition 10, we 
have 

{U - Afa) = E X (-Af ) = 4cisin 2 j^j . 

Furthermore due to the convexity of G(t) and Jensen inequality we can estimate for p G 9Jti(Aj) 

Y,v{x)S(p{x),t) < G(t)/t-£ d GU- d tJ2p^))/t= S(r d ,t), 

xgA \ xeAf I 
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respectively 



□ 



Proof of Proposition 8. Combining the Feynman-Kac representation of {u(t, 0)} in (2) and N(t) 
in (6) (u(t,0) > N(t) is obvious. With the hitting time defined by t\(x s ) := inf s >o [x s G A c ] . 
and averaging with respect to A due to ergodicity, we obtain the lower bound 

V u (x s )ds ) 5 x (x t ) l TA>t ) 



N(t)> \A\~ l J2(w\exp(- [ 
, c a L V Jo 



Choosing I adequate and applying Proposition 11 we obtain 



N(t) > exp lG(t)-tM X 7(t)0- + o(l)) 



3. The upper bound on a box 



□ 



In this Section we assume that the hypotheses of Theorem 2 are satisfied. We want to prove 
the following upper bound. 

Proposition 12. Denote by Ei(H^) the principal eigenvalue of the Hamiltonian on 
A = Aj(0), / G N defined in (4)- With the Faber-Krahn constant cfk, 7 = cfk /(12vr) 2 , 



and 



xt(t)-= max min 2d(l - 2Vl - h),-f/2 + (1 - h)S((l - h),t)) 

0.5<h<l L 

Xt(t) = 7sm 2 (f ^) + S((A£)- d ,t)/4, 



£ > 2, we have 



(expi-t Ei(Hg))) < exp (c(t) -trdxf(t) + C\A\\ . 



Trying to transfer the proof of Proposition 8 to Proposition 12 we observe two difficulties. We 
have to control each p G 9Jii(Z d ) and it is not possible to interchange in (33) the expectation 
value and the supremum. To solve the first problem we restrict to the cube A = A/(0) := 
{x G Z d : \x\oo < 1} with I G N to be choosen in the next section and define the following 
classification of p G SDTi (A) . 

Definition 13 (Classification of the occupation time measures ). 
With 7 as in Proposition 12 we define 

5(1) := [p G arti(A) : 7 sin 2 (~) < (^| - A^) < 2d)}, 

and, if £ > 2, 

m ■= {p G 3(1) C ■■ 7sm 2 (f ^i) < (VP\ ~ *aVp) < 7^ ( J)} • 
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The classification of p £ SJTti (A) in Definition 13 is reminiscent of the contraction principle in 
large deviation theory [12]. The energy form of the Laplacian plays the role of the projection 
map while $(£) correponds to the resulting value set. 

We have to deal the case I = 1 separately to prove the asymptotics for the classical regime, i.e. 

(u(t, 0) = exp(G(t) - 2dt{\ + o(l))). 
Moreover the argument below is prototypical for general I. 
Definition 14. Fix 1 > h > 0.5. We define 

5^(1) := {p £ 5(1) : max p(x) > h}, 

r h (l) :=$(l)\3h0-)- 
Let us first deal the situation when p £ 5(1) is "almost" a 5-peak. 
Lemma 15. 

( sup exp (t[(y/p\ Afty/p) + (y/p\ Vuy/p)])) < |A| exp (G(t) - 2dt(l - 2\/Y^h)) . 
Proof. Define x\ by p{x\) = max x ^\p(x). Then 

(y/p\ -A%y/p)> (Vl>fa)- VpV)) 2 > 2d(l-2y/T=h), 

\x\— x\=l 

respectively 

( sup exp(-t[{y/p\A^y/p) + {y/p\V uy /p)])) < |A| exp ( G(t) - 2dt{\ - 2x/T=h)) . 
pe? h (i) v ' 



□ 



Lemma 16. 



( sup exp(-t [(VP | Aj>y/p) + (jp\V uy /p)])) 
pe^(i) 

< exp (G(t) - 7t/2 - (1 - fr)iS((l - /i), t))) 



Proof. Define xi, xi by 



V w (xi) = minVuix), 

V U) (x2) = min V u (x). 
xeA\{xi} 
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Then 

( sup exp ( -^p(x)V u (x)t I) 

< X] ( sup exp ~ ^Pi^t (Vu( x ) - V u (yi)) : yi = x x ) 

< V] (sup exp(-V u (yi)t - (1 - p{yi))t {V UJ (y 2 ) - V u (yi))) : = x u i = 1,2) 
yieA pess;(i) 

y 2 GA\{j/i} 

< (exp - (1 - ^^,(2/2))) 

2/1 eA 
y2eA\{i/i} 

< |A| 2 exp(G(/it) + G((1 - fc)t)) 

< |A| 2 exp(G(t) - (1 - /i)t5((l - /*),*)), 

respectively 

( sup exp(-t [(VPl A A 5 VP) + (v^l^v^)])) 
peS c h (i) 

< exp {--ft/2) ( sup exp (-t(y/p\ V uy /p ))) 

pe^(i) 

< | A| 2 exp (G(t) - "ft/2 - (1 - h)tS((l - h),t))) . 



□ 



Combining Lemma 15 and 16 we have 

Corollary 17. With x± (~t) defined in Proposition 12 we have 

( sup exp(-t[(^p\A%^p) + (^p\Vu^p)})) < \A\ 2 exp(G(t)-txt(t)). 

PG3(1) 

We now discuss £ > 2. 

Definition 18. Let h > and p € 9Jti(A). The level set is denoted by 

U h (p) :={xeA: p(x) > h} . 

For £ > 2, n G No we define 

^,n) := {pe5(l) = n} . 

The first step to prove an analogue of Corollary 17 for $(£, n) is the following generalization of 

y p(x)W(x) > mmW(x) } p(x) = mmW(x). 
^— f xeA f xeA 

xeA x-GA 
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Lemma 19. Let {W(xj)}j = i t iju be an arrangement ofW : A — > [0, oo) according to the size, 
that is 

W(xi) = minW(x), 

W(xj) = min W(x), 

x£A\{xi,..Xj_i} 

p £ $(£,n), and suppose that 

L^ _1 |bXc/^(p)IIJ : = max{m £ N : m < hT 1 III- 

Then, 



Y p(x)W(x) > < 



\\PXu-(p)\\i < h, 

"-z^j=i vv{x n+ j) otherwise. 



Proof. We only discuss ||pxt/=(p)l|i — h. With {xj} as denned above, ||xc c (p)Pl|l — h\U£(p)\ 
and J = n + j we have 

i h ~ 1 \\Xuf i (p)P\\l\ i h ~ 1 \\Xuc(p)P\\l\ \U£(p)\ 

Y h - \\xvfa)p\\i = Y p ( xj } + p ( xj }> 

7=1 7=1 i=l h ~ 1 \\xuf i (p)P\\i\+ 1 



respectively 



[h ^{xu^phl \U£(p)\ 

Y (h-p(xj)) < Y P( X J)- 

i =1 i=\h~ 1 \\xim P )Ph\+ 1 



As a consequence of W : A — > [0, oo) we have 



l h 1 IIXC7^(p)PlllJ 

Y W(xj)(h-p(xj)) 

7=1 

L' l » -1 llxi7=(p)PlliJ 

< W(x n+lh - H phi ) Y (h-p(xj)) 

7=1 

<w(x n+lh - HxnMPhi ) Y p( x ^ 

j=ih- 1 \\xuf t (p)P\\l\+ 1 

M(p)\ 

< Y p{xj)W{xj). 

j=\ h ~ 1 \\Xuf i ( P )P\\i\+ l 
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Thus, 

E v{x)W{x) 

l h ~ 1 \\Xuc(p)P\\l\ 

= h e w{ ~ x j) 

i=i 

^"'IIx^wpIUJ \u£(p)\ 

E W(xj)(h-p(xj))+ J2 p(xj)W(xj) 

L ft_1 !lxc7C(p)PlliJ 

> h e 

□ 

The second ingredient to prove an analogue of Corollary 17 is a lower bound of ||xt/^(p)Plli f° r 
p £ n) proven in Corollary 22. To attain this goal some intermediate steps are necessary. 

Lemma 20. Suppose p G 9Jti(A), h > 0. Define pj L : Z d -»■ E &y 

f(\^M-^) 2 ifxeU h {p), 

Phyx) := < 

I otherwise. 



(VPl-Afvp) > (v^l - Ag j(p) V^ 



Then 



Proof. As a consequence of the definition of 50? i (A) in (34) we have 

(Vp|-a£vp) 

= \ E (VpM- v^)) 2 + E E (Vp@)- vpM) 2 

x,yeA, x6A yeA c , 

|x— i/|=l |as— J/|=l 

> \ E (vpM-v^)) 2 + E E (V^)-^ 

x,yeU h (p), xeU h ( P ) y€U£( P ), 

jas— w|=l |a;-2/|=l 

> ^ E (v / ^M-v / ^)) 2 + E E (V^)-^) 2 

x,y£U h (p), xeU h (p) yeU^ip), 

\x-y\=l \x-y\=l 

D 



Ph\ ~ &U h ( P )VPh 



□ 
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Proposition 21 (Faber-Krahn-inequality[5]). Suppose U C Z d , \U\ < oo. Then there is a 
constant < cfk < 2d, s.t. the principal eigenvalue Ei(-A^) of the discrete Dirichlet- 
Laplacian — A® is bounded below by 

Ei(-Ag) > c FK \U\- 2 ' d . 

As a consequence of the Faber-Krahn-inequality we find a lower bound of the probability mass 
outside the level set. 

Corollary 22. Suppose that £>2,p£ $(£,n) and 7 = cfk/(12tt) 2 . Then 

Proof. By Definition 13, Lemma 20 and Proposition 21 we have 

7 sin 2 (J) > (VPl-Afvp) 

* E ^- A U h(p) )\\V^\\ 2 2 
> c FK n- 2ld \\^\\l, 



respectively 



and 



fell, <^n 2 l d sin 2 (l)<^-ry/- 



|Xt/j,(p)Pl|l < l|XDfc(p)[^- V^+V^JIll 

< \\\XU h (v)[Vp-^\h + \\XU h (p)^\\\ 



2 

2 



respectively 

\\Xu^ P )P\\i = 1 - \\Xu h {p)P\\l > 



In the next Lemma we choose the height h in dependance on i. 

Lemma 23. Suppose that £>2,h = {U)~ d and p G $(£, n), n < h' 1 = (A£) d . Then 

2<N r .= (U) d /2 - 1 < n+ L/» _1 ||x^(p)PlliJ, 

and 

h(N e - 1) > 1/4. 



□ 
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Proof. First observe 

h(N e - 1) = 1/2 - 2h > 1/4. 

and nh = 4~ d n£~ d < 1, respectively n l l d l~ l < 4. Combining the assumptions of Lemma 23, 
Corollary 22 and 7 = cfk /(12vr) 2 denned in Proposition 12 we get 

n/i + /i[/i _1 ||x^(p)P||il > nh+ \\xU£(p)P\\i 

( / 2 \ 1/2 \ 2 

>l_((_T^j n Vd r i + (h n flA +nh 

>l-^n^-2(^-) 1/2 n^ 
Ac FK \4:C FK J 

2 / 2 \ 1/2 

7"7T / 77T N 



respectively 



> 1 - - — 4-4 

cfk V cf_^ 

> 1/2, 



n+ L^ _1 ||x^(p)PlllJ > ^ > 2. 



□ 



The next proposition generalizes the argument given in the proof of Lemma 16. 
Proposition 24. Suppose £ > 2 and as defined in Lemma 23. Then 

( sup exp (-tJ2 p(x)K,(x))) < |A|( |A [ ) exp (c(t) - 4~ 1 ^((4^, t) 

Proof. Assume {V^*}^^..,^! is an arrangement of V : A — > R according to the size, that is 

V* = V u (xi) = mmV UJ (x), 

x£A 

Vj = V u (xj) = mm V u (x). 

xGA\{xi,..Xj-i} 

and {p*}i=i |A| is defined by 

Pi =p{Vi) = maxp(y), 
Pi=p{Vi)= max p(y). 
An elementary induction argument gives 

|A| 

xgA j=i 

respectively 

( sup exp -tV p(x)K,(x) ) < ( sup exp -fV ). 
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Suppose h = (4l)~ d and N t > 2 as in Lemma 23. Observing V* - V* >0,p*> h, j = 1, ..,n, 
Lemma 19 and Lemma 23 give for p £ n) 

|A| n |A| 

Y /P *v* = v 1 * + J2p*A v /- v i)+ E p*A v ?- v i*) 

3=1 j=l j=n+l 

n +\ h ~ 1 \\xuf i ( P )Ph\ 
> V* + h i V j ~ y i) 

3=1 

>V{ + h ^2(V* - V?) 

3=1 

= (l-h(N e -l))v; + hJ2 V*, 

3=2 

respectively with Lemma 23 



( sup exp j p(x)V w (x) ) 



< (expHl-^-l))^-^^^*)). 



3=2 

= E E (eM-(l-HNi-l))V w (x)-hY, V u (y)) :x = xl,W = {x* 2 ,..,x* Ne }). 

z€A wcA, y£W 
\W\=N f -l 

< E E < ^(-c 1 - w - i))^^) - ^ E ^(^))>- 

xeA wcA, y£W 
\W\=N e -l 

= l A l _ i) < exp (-^ - W - 1) W(0))> II <«p(- W(0))>. 
Integration with respect to the single site potential and Lemma 23 give 



( sup exp -*E 

p6j(*,») V x-GA / 



- |A| Cvi - i) exp " h{Ne " m + {Ne ~ l)G{th)) 

- |A| (ivi - i) exp " h(JV * " 1)ts{h > t]) 

- \ A \(J^_^w(G(t)-tS(h,t)/4,). 



□ 



Proof of Proposition 12. As a consequence of the Faber-Krahn-inequality exists I = cl s.t. 
£0li (A) C U{=i Choosing h = {Al)~ d as in Proposition 24 we observe |Ufc(p)| < {U) d < 
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(Al) d . Corollary 17 gives 

<exp(-i E 1 (H»))) 

< ( sup exp(-t[(vp| -A^VP) + (v^|K,^)])) 

pej(i) 

i {U) d 
£=2 n=0 PeW") 

< \A\ 2 exp (G(t)-t X f(t)) 

+ g exp (- 7 sin* (__jj g (^«P<-* ( VPI V.^)». 

As a consequence of Proposition 24 we obtain 
<exp(-ti?i (#£,,))> 

< |A| 2 exp (G(t) - t X t(t)) +C\A\J2 (J M x ) (<?(*) - ttf(f)) 

< |A| 2 exp(G(t)-t X +(t))+|A| 2 exp (^) - t mf X +(t)) £ ^ |A| 
<ex P (G(t)-tMxt(t) + C\A\y 

4. The proof of Theorem 2 



Nt-1 



□ 



To prove Theorem 2 we have to choose the side length I = l(t) of the cube A = A;(0) s.t. the 
error in Proposition 12 and the error resulting from boundary conditions are negligible, that is 



c^tf + c 2 ti(ty 2 o (tM X t(t)) 



To solve this problem we have to compute inf^g^ x/(i) in dependance on the single site distri- 
bution, respectively the corresponding cumulant generating function. In Theorem 2 we assume 
G(t) < oo, t > and G(t)/t E II 3 with auxilary function g E R p , p £ [—1, oo) and g-index c g , 
i.e. we have 

S(X, t) = m_m_t-±°° Cghp {X)g{t){l + o(l)), 
A E (0, 1]. As discussed in [3], p. 128, h p (X), X E (0, 1] has the representation 



MA)= / u^du + o(l) = \; log{X } / P = (40) 
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i.e. to estimate inf^g^x/W we have to distinguish three cases. Let us first discuss single site 
distributions exhibiting Lifshitz tail behaviour. 

Lemma 25. Suppose G(t) < oo, t > and G(t)/t £ H g with auxilary function g £ R p , 
p £ [—1,0) and tg(t) — > oo in the limit t — > oo. Denote the optimal length by 



t(t) := g(t) l ^ d P~ 2 \ 



Then 



C 1 t(t)- 2 < M X }{t) < inf X J(t) < C 2 t(t)-\ 



G\, C2 > 0, t sufficiently large. 

Proof. To prove the lower bound compute the corresponding continuous minimization problem. 
As a consequence of £*(t) — > oo the upper bound is given by 

inf Xt - (t ) < W (Ijt^ri) +s«rw<.t) £ WW" 1 - 

□ 

In the next lemma we discuss p = defining the borderline between the classical and the quan- 
tum regime. It contains the single peak case (£*(t) = 1), the double exponential distribution 
(£*(t) ~ and the almost bounded single site distributions (£*(t) — > oo, t — > oo). 

Lemma 26. Suppose G(t) < oo, t > and G(t)/t £ Tl g with auxilary function g £ Rq. Then 

'2d c g g(t) > 2vr 2 , 

8dc g g(t) + 2dc g g{t) log (2ir 2 /{c g g{t))) c g g(t) < 2ir 2 



(41) 



and 



inf xt(t) > 

££N e V ' 



'2d-4d(c g g(t)) 
min 



-1/2 



/9 rfc g g(t) rfc g g(t) , / 64 7 7r 2 \ 

Proof. Approximating sin(x) and choosing 

t(t) = max 



c g g{t) > 2e 2d + 7vr 2 /2d, 
c g g(t) < 2e 2d + -fir 2 /2d. 



(42) 



i 

\c g g{t)J 

we obtain (41). The starting point to prove (42) is 

^7,-2 , d 



inf Xi (*) > min Xi inf 



£e[2,oo) 



+ -c ff5 (t)log(4£) 



Computing the infinum with 



t(t) = max 



/ 4vr 2 7 \ 1/2 ' 
V 5^(0 J 



(43) 



(44) 



the lower bound (42) is then a consequence of some elementary, but painful calculations. □ 
Finally let us discuss the one-peak case 
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Lemma 27. Suppose G(t) < oo, t > 0, G(t)/t G II 3 mi/i auxilary function g G i? p , p G (0, oo). 
T/ien f (i) = 1 and 



2d = inf x7(t) > inf vt(t) > 2d 1 - 2o(t) 



-1/2 



(45) 



t— >oo 



for t sufficiently large. 

Proof. By assumption we have 

S(X,t) =c gP - 1 {l-\P)g{t) + o{l) -> oo. 

Choosing h = l- l/g(t) > 1/2 we obtain (45). □ 

Corollary 28. Suppose G(t) < oo, t > and G(t)/t G II 5 wii/i auxilary function g(t) = 
t p 9o(t) ^ -Rp> P ^ [ — 1)°°)j 50 £ -Ro anc ^ g-index c g . Furthermore if p = — 1 we assume that 
9o(t) *-^>° oo. Defining I = l(t) := |~a(i)£*(i)] wi/i £*(t) as Lemma 25 -27 and 

' 'go(t)v&w p=-l, 
i+p 

t 3(^2) p g (-1,0), 

p G [o, oo) 



a(i) 



+ K 2 tr 2 = o (t inf j • 

\ em 1 J 



(46) 



Proof. As a consequence of inf^g^x/W > we can estimate 

K x l d + K 2 tr 2 



< 



< 



K x l d + K 2 t/" 2 



< Ka{t) d t- l t{t) d+2 + Ka{t)~ 



The case p G [0,oo) is obvious. Suppose p G [-1,0). With a(t) ->■ oo and f (t) = #(i) 1//(<ip ~ 2) 
we have 



< 



< Ka(t) d t d p-2 5f (i) dp " 2 + Ka(ty 2 ->■ 0. 



d+2 



□ 



The next lemma is a slight modification of Proposition 4.4 in [4] to deal the bounded and the 
unbounded setting simultaneously. 

Lemma 29. Suppose L := L(t) = [~ilog(i)~|. Then 

(exp (-tEi < exp f G(t) - t inf x +(t)(l + o(l))) . 

Proof. Lemma 4.6 in [4] says, there is a constant C > such that for every I G N, there is a 
function <pi : Z, d — > [0, oo) with the following properties: 

(i) (pi is /-periodic in every component, 

(ii) INIoc < C/l 2 , 

(iii) For any potential V : Z rf — > [0, oo) and any L > 1/2 we can estimate 



Ei((-A + ^ + &)aJ > min E 1 l(-A + V)° ix) 
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We define 



V u :Z d ^[0,oo), 



V u (x) 



- mm xe \ 2L V w (x) x 6 A 2 l, 



otherwise. 
With H u := — A + V w + (pi we obtain as a consequence of the result above 
E x (fffj = inf - + min K,(x) 



/ei 2 (A_L) L 
lb=l 



>£i < - CP 2 + min V u (x) 



x€A 2L 



respectively 



(exp(-^i (#£J)> 



< exp(CT 2 t)(exp ( (HP) - t min K,(x 

/ xgA 2 z, 



< exp(CT 2 t)( max exp f-i ^ \ H^ x) 



< exp(Cr 2 t) ^ (exp (-t E x (h 



D 

A;(z) 



< 2 d \A L \ exp(Cr 2 t)(exp (-i ^ (ijf ;(0) 

Choosing / = = [a(t)£*(t)] as in Corollary 28 we have l(t)/2 < L(t) for t sufficiently large. 
An application of Proposition 12 gives 

(exp (-t Ex < 3 d \A L \ exp \G(t) - tiri X t(t) + Cxl(t) d + C 2 tl(t)~ 2 ^j . 

Lemma 29 is now a consequence of Corollary 28. □ 
Proof of the upper bound of Theorem 2: Choose L = \t log(t)] . Then 

(n(t,0)) = (E° r 

= <E° ex p("^ V„{x.)te\ 1 

The second term can be estimated by 

<E° exp (j V u (x s )ds^\ 1 

Applying the estimate of the hitting probability 

IP [r L (x t ) < t] < 2 d+1 exp (-L (log (L/(td)) - 1)) 



■ / V u {x s )ds 



TA L <t 



+ (E° 



exp f- J V u) (x a )dsj ^TA L >t 



< |A|(exp(-tE 1 «J)). 
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Lemma 2.5) we have 



exp (J V w (x s )ds\ 1 



<E° i exp(G( 

< 2 d+l exp (G(t) - t log(t) (log (log(t)/d) - 1)) 
Finally as a consequence of Lemma 29 we have 

(u(i,0)> < exp(G(t)-tnif Y+(t)(l + o(l)) 



in the limit t — > oo. □ 
Proof of Corollary 4 ■ Combining Theorem 2 and Lemma 25 we obtain 

— 2 —2 

-C 1 tg(t)'^ < logiV(t) < log(it(t,0)) < -C 2 tg(t)^2. 

Applying Corollary 31, i.e. the limit of oscillation version of de Bruijn's Tauberian theorem in 
Appendix 2, we obtain upper and lower bounds of the IDS in the limit E \ 

Cimf[Et-tg(t)^](l + o(l)) < log N(E) < C 2 inf [Et - tg(t)^](l + o(l)). (47) 



Observing that 



([3], P-44) the minimizing time t* of the Legendre transform satisfies with C > 1 



0(f) = (t*)Pg (t*) CEP-tpy* 



(48) 



(49) 



As a consequence of (13) we can apply the inversion formula for regularly varying functions 
stated at the end of Appendix 1 and obtain 



f CE^ dp y 2p g (cE {2 ~ dp ^ 2p 



respectively 
inf 



Et - tg(t)*p 



-i/p 



go 



-(C-l^-i+^^o (CE^-^I 2 P 



-l/p 



□ 



Proof of Corollary 5 -7 . Combining Theorem 2, Lemma 26 and Lemma 27 gives Corollary 5, 
i.e. the estimate 



G(t) - 2dt X *_(t)(l + o(l)) < log iV(t) <log(n(t,0)) < G(t) - 2dt X *+(t)(l +o(l)) (50) 



with 



fl c 95 (t)>2vr 2 , 
\4c s5 (t) + c 9 5(*) log (27T 2 /(c 9 5(0)) <*s(t) < 2^ 2 
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and 

■ 1 - 2(c g g(t)r 1 / 2 c g g(t) > 2e 2d + 7 7r 2 /2d, 



mm 



7/(4d), ^# + ^fM log ( ) ^(t) < 2e M + 7^/2d. 



Suppose now p > 0. Then we have — > oo and 1 = x*-{t) > X+{t) = 1 — i n the 
limit t — > oo. Applying Corollary 33 in Appendix 2, that is the limit of oscillation version of 
Kasahara's Tauberian theorem, we obtain in the limit E — > — oo 

log N(E) = M[{E - 2d)t + G(t)](l + o(l)) = -/(£ - 2d)(l + o(l)). (51) 

In the double exponential setting we have c g g(t) — > c g and X-W as wen as X+(t) converge to 
constants. Corollary 6 follows now from the analogue of Corollary 33 in the double exponential 
case proven in [24] (see also Appendix 2). □ 

Appendix 1: Regular varying functions 

Regularly varying functions as introduced in Definition 1 are a generalization of g(t) = t p . 
Their characteristic trait is 

g(Xt)/g(t) t ^°X p (l + o(l)) 
for all A > 0. Sometimes it is convenient to transfer attention from infinity to the origin. Thus 

if g > o, 

g(XE)/g(E) E ^° A'(l + o(l)), 

we say g is regularly varying at the origin with index p, g € R p (0+). This is equivalent to 
g(l/E) £ R- p , [3], p. 18. If p = 0, then g is said to be slowly varying. Regularly varying 
functions are a generalization of g(t) = t p in the sense that g S R p implies g(t) = t p go(t) with 
go G Ro [3], Theorem 1.4.1. 

One problem is the inversion of regularly varying functions. Theorem 1.5.12 in [3] says, if 
g € R p with p > 0, then exists an asymptotic inverse g~ l € R\/ p with 

Up to asymptotic equivalence g -1 is unique. A corresponding result in the case p < with 
g~ l £ i?i/p(0+) can be deduced from Theorem 1.5.12. To obtain an explicit expression for 

it 

the asymptotic inverse, we introduce the de Bruijn conjungate g(f ([3], p. 29), i.e the slowly 
varying functions go G i?o satisfying 

2o(t)5o # (^oW) 1, 9o # (t)5o(^ # (i)) ^° I- 
Again, up to asymptotic equivalence g$ is unique. Suppose now that 

5o(^o(t) 1/p )/5o(t)'^°l (52) 

holds for some p ^ 0. As a consequence of Corollary 2.3.4 in [3] we have (g^ p )^ ~ g and 
if E ~ t p 5o(i) then t ~ E 1 / p g (E 1 / p )~ 1 /p in the limit t, £ -> oo if p > 0, respectively t ->■ oo, 
E 1 \ 0, p < 0. An example satisfying (52) is given by go(t) = log(t). The interested reader is 
encouraged to control the statements above for g(t) = t p log(t). 
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Appendix 2: Tauber theory for Laplace transforms 

In Appendix 2 we collect some results about Tauberian theorems. Given the asymptotic 
behaviour of the Laplace transform 

N(t) := J e" A * dv{\) < oo (53) 

in the limit t — > oo, the problem is to reconstruct the behaviour of the distribution function 
v(E) = v[Eq,E] in the limit E \ Eq = sup{S : v{E) = 0}. This is a very common problem 
in statistical physics, respectively probability theory, and the aim is a characterization of the 
asymptotics in terms of the Legendre transform. 

log (v(E)) E ^ Eo inf [Et + log(iV(t))] (54) 

(e.g. [26], Thm. 9.7). If the exponent of the Laplace transfrom is a regularly varying function, 
very explicit statements are possible. In the bounded setting one has de Bruijn's Tauberian 
theorem ([3], Thm. 4.12.9). The corresponding result in the unbounded setting is Kasahara's 
Tauberian theorem ([3], Thm. 4.12.7). 

To apply de Bruijn's, respectively Kasahara's Tauberian theorem one has to know the exact 
asymptotics of N(t). Sometimes this is a practical problem, because only upper and lower 
bounds of the Laplace transform are available. But if these bounds are in good agreement, 
we can apply the so called limit of oscillation theorems. We start by discussing the limit-of- 
oscillation version of de Bruijn's Tauberian theorem ([2], Thm. 0) 

Theorem 30. Let v be a measure on (0, oo) whose Laplace transform N(t) satifies (53). By 
v(E) we denote the distribution function of v. If a < —1, tjj G R a (0+), put '■= At/^A) G 
R a+ i(0+). Suppose 

-Bx< lim inf A log N(ip(X)) < lim sup A log N(tp(X}) < -B 2 (55) 

*\0 A\0 



for B 1 ,B 2 >0. Then 



-Ci < lim inf Au(l/<£(A)) < lim sup Au(l/0(A)) < -C 2 (56) 
*\o A\0 



with Ci, C 2 > 0. 



Corollary 31. Let v be a measure on (0,oo) whose Laplace transform N(t) satifies (53) and 
denote by v(E) its distribution function. Suppose f G R p , < p < 1 and 

-B 1 f(t)(l + o(l))<log(N(t)))<-B 2 f(t)(l + o(l)) (t^oo) (57) 

for Bi,B 2 >Q. Then 

d inf [Et - /(i)](l + o(l)) < log (v(E)) < C 2 inf [Et - /(i)](l + o(l)) (E \ 0). (58) 
with Ci,C 2 > 0. 

Proof. Again we encourage the interested reader to control the statements below for the special 
case f(t) = t p , p€ (0, 1). With ip~ l {t) = l/f(t) G R^ p inequality (57) becomes 

-Bi < liminf tp' 1 (t) log N (t) < limsup ^~ l {t) log N(t) < -B 2 , 
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respectively with A = 1/t and ^>(A) G R-i/ p (0+) 

-B x < limsup A log iV(V>(A)) < liminf A log iV(^(A)) < -B 2 . 

A\0 A \° 

As a consequence of Theorem 30 and (f>(\) = A^(A) G Rp-i (0+) we have 

p 

-Ci < liminf \v(l/4>(\)) < limsup Av(l/0(A)) < -C 2 . (59) 
A \o a\o 

Setting = 1/^(A) we obtain 

-C 1 /^- 1 (l/£;)(l + (l)) < log (v(E)) < -C 2 /<P'\l/E)(l + o(l)). (60) 

Let us now prove (58). Without restriction we can choose a differentiable version of f(t) [3]. 
Starting with 

f 9 1 



E 0- 1 (p/E) 
we have 

p/*~*(w)-*(w)w- (61) 

Inversion of ?/> gives ^ _1 (i*) = p/(Et*), respectively 

= £ - iW) = £ " ^C*)/'* E ~° £ " >> ■ < 62 > 

We obtain 

- /(t)1 = £ IFmW - ^ - ^"^-'(l/S). (63) 

Combining (60) and (63) we end up with 

^— injp-/(t)](l + (l)) < log (v(E)) < 2P i _ p mf[Et-f(t)](l + o(l)). 

□ 

The corresponding result in the unbounded setting is the limit-of-oscillation version of Kasa- 
hara's Tauberian theorem, ([13], Thm. l(ii)). 

Theorem 32. Let v be a measure on K whose Laplace transform N(t) satifies (53) and denote 
by v(E) its distribution function. If < a < 1, ip G R a , put <f>(t) = t/ip{t) G R\- a - Suppose 

Bi < liminf t~ l log N(ij)(t)) < limsup t' 1 log N(tp(t)) < B 2 

for B t ,B 2 >0. Then 

d < liminf E~ l logv(-(p(—E)) < limsup E" 1 log v(-<p(-E)) < C 2 

E^-oo E->-oo 

with Ci,C 2 > 0. 
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Corollary 33. Let v be a measure on R whose Laplace transform N{t) satifies (53). Denote 
by v(E) the distribution function of v. Suppose f € R p , p > 1 and 

B 1 f(t)(l + o(l))<log(N(t)))<B 2 f(t)(l + o(l)) (i-Kjo) (64) 

for B 1 ,B 2 >0. Then 

d inf [Et + f(t)](l + o(l)) < log (77(E)) < C7 2 inf [Et + f(t)](l + o(l)) (E -> -oo). (65) 

Proof. Corollary 33 may be proved by a method closely analogous to that used in Corollary 
31. □ 

Remark 34. 

(i) The constants C\,C% > in the results above are explicetly computable. As a con- 
sequence if B\ = B 2 then C\ = C 2 , i.e. we obtain de Bruijn's (['!], Thm. 4.12.9), 
respectively Kasahara's Tauberian theorem ([3], Thm. 4.12.7). 

(ii) The key idea linking Laplace and Legendre transform is the concept of the relevant 
energy interval. In the asymptotic limit for every time t exists an energy E = E(t), 
s.t. the behaviour of the Laplace transform is determined by a small energy interval 
around E. Using the Cramer-transform ([12], p. 7) the idea above is used in [4] to prove 
de Bruijn's Tauberian theorem, repectively in [24] not aware of [13] and [2] to prove 
the corresponding limit of oscillation versions. 

(hi) While in general no Tauber theorem exists if N(t) £ R\ it is possible to transfer 
the limit of oscillation argument via Cramer-transform discussed in (ii) to the double 
exponential regime [24]. Suppose B\,B 2 > 0, f(t) = c g t\og{c g t) — c g t and 

B 1 f(t)<logN(t)<B 2 f(t). 
Then exists C\.C 2 > s.t 

Ci inf [Et + /(*)] (1 + o(l)) < \ogv(E) < C 2 inf [Et + f{t)] (1 + o(l)). (66) 
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